The Dark Matter Annihilation Signal from Dwarf Galaxies and Subhalos 
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Dark Matter annihilation holds great potential for directly probing the dumpiness of the Galactic 
halo that is one of the key predictions of the Cold Dark Matter paradigm of hierarchical structure 
formation. Here we review the 7-ray signal arising from dark matter annihilation in the centers 
of Galactic subhalos. We consider both known Galactic dwarf satellite galaxies and dark clumps 
without a stellar component as potential sources. Utilizing the Via Lactea II numerical simulation, 
we estimate fluxes for 18 Galactic dwarf spheroidals with published central densities. The most 
promising source is Segue 1, followed by Ursa Major II, Ursa Minor, Draco, and Carina. We show 
that if any of the known Galactic satellites can be detected, then at least ten times more subhalos 
should be visible, with a significant fraction of them being dark clumps. 



INTRODUCTION 

A decade has gone by since the emergence of the "Miss- 
ing Satellite Problem" [l|, Q , which refers to the apparent 
discrepancy between the observed number of Milky Way 
satellite galaxies, 23 by latest count d 0, H, S, 0, i, B [13, 
11, l2|, and the predicted number of dark matter (DM) 
subhalos that should be orbiting in the Milky Way's halo. 
The latest cosmological numerical simulations [ij, [13, [ill 
resolve close to 100,000 individual self-bound clumps of 
DM within the Galactic virial volume - remnants of the 
hierarchical build-up of the Milky Way's DM halo. A con- 
sensus seems to be emerging that this discrepancy is not 
a short-coming of the otherwise tremendously successful 
Cold Dark Matter (CDM) hypothesis [il[i3|, but instead 
reflects the complicated baryonic physics that determines 
which subhalos are able to host a luminous stellar compo- 



effects on strong gravitational lensing (see Zackrisson & 
Riehm's contribution in this special edition), or from the 
kinematics of stars orbiting in DM-dominated potentials 
3l|, or from perturbations of cold stellar structures like 



globular cluster tidal streams 32, 3 ^ 34^ 351 o r the heat 
ing of the Milky Way's stellar diskjHllillii] , although 
all of these are also worthwhile approaches to take. 

The only trouble with the DM annihilation signal is 
that so far there have been no undisputed claims of its 
detection. Recently there have been several reports of 
"anomalous" features in the local cosmic ray flux: the 
PAMELA satellite reported an increasing positron frac- 
tion at energies between 10 and 100 GeV [s^, where 
standard models of cosmic ray propagation predict a de- 
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creasing fraction; the ATIC [40| and PPB-BETS 
balloon-borne experiments reported a surprisingly large 
total electron and positron flux at ~ 500 GeV, although 
nent and which aren't [H [l9|, ll^, IMI, [H III] . If this recent Fermi ^ and H.E.S.S. data ^ appear to be 



explanation is correct, then an immediate consequence 
is that the Milky Way dark halo should be filled with 
clumps on all scales down to the CDM free-streaming 
scale of 10-12 iq-^Mr (li,!!!. At the moment there 



is little empirical evidence for or against this prediction, 
and this has motivated searches for new signals that could 
provide tests of this hypothesis, and ultimately help to 
constrain the nature of the DM particle. 

One of the most promising such signals is DM an- 
nihilation 1271]. In regions of sufficiently high density, 
for example in the centers of Galactic subhalos, the DM 
pair annihilation rate might become large enough to al- 
low for a detection of neutrinos, energetic electrons and 
positrons, or 7-ray photons, which are the by-products 
of the annihilation process. This is one of the few ways 
in which the dark sector can be coupled to ordinary mat- 
ter and radiation amenable to astronomical observation. 
Belying its commonly used name of "indirect detection" , 
DM annihilation is really the only way we can hope to 
directly probe the dumpiness of the Galactic DM dis- 
tribution. One could argue that it is a more "direct" 
method than trying to constrain DM dumpiness from its 



inconsistent with it. Either of these cosmic ray anoma- 
lies might be the long sought after signature of local 
DM annihilation. However, since the currently avail- 
able data can equally well be explained by conventional 
astrophysical sources (e.g. nearby pulsars or supernova 
remnants) , they hardly provide incontrovertible evidence 
for DM annihilation. The next few years hold great po- 
tential for progress, since the recently launched Fermi 
Gamma-ray Space Telescope will conduct a blind sur- 
vey of the 7-ray sky at unprecedented sensitivity, energy 
extent, and angular resolution. At the same time. At- 
mospheric Cerenkov Telescopes, such as H.E.S.S., VERI- 
TAS, MAGIC, and STACEE, are greatly increasing their 
sensitivity, and have only recently begun to search for a 
DM annihilation signal from the centers of nearby dwarf 
satellite galaxies [13, [H, [H, [13, [H ■ 

The purpose of this paper is to provide an overview 
of the potential DM annihilation signal from individual 
Galactic DM subhalos, either as dwarf satellite galaxies 
or as dark clumps. It does not cover a number of very 
interesting and closely related topics, which are actively 
being researched and deserve to be examined in equal 



2 



detail. These include the diffuse flux from Galactic sub- 
structure and its anisotropics [e.g. 4^, 5^, 51 1, the relative 
strength of the signal from individual subhalos compared 
with that from the Galactic Center or an annulus around 
, the effect of a nearby DM subhalo on the am- 
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plitude and spectrum of the local high energy electron proximately 



Hubble constant at cosmological scale factor a) exceeds 
the annihilation rate T{a) = n{av), i.e. when DM par- 
ticles can no longer find each other to annihilate. The 
co-moving number density of DM particles is then fixed 
at a "freeze-out" temperature that turns out to be ap- 



and positron flux 53 , |54| | , and the role of the Sommer- 



feld enhancement [55| on the DM annihilation rate and 
its implications for substructure signals 

This paper is organized as follows: we first review the 
basic physics of DM annihilation, briefly touching on the 
relic density calculation, the "WIMP miracle" , DM par- 
ticle candidates, and, in more detail, the sources of 7-rays 
from DM annihilation. In the following section we review 
what numerical simulations have revealed about the ba- 
sic properties of DM subhalos that are relevant for the 
annihilation signal. We go on to consider known Milky 
Way dwarf spheroidal galaxies as sources, using the Via 
Lactea II simulation to infer the most likely annihila- 
tion fluxes from published values of the dwarfs' central 
masses. Next we discuss the possibility of a DM annihi- 
lation signal from dark clumps, halos that have too low 
a mass to host a luminous stellar component. Lastly, we 
briefly discuss the role of the substructure boost factor 
for the detectability of individual DM subhalos. 



DARK MATTER ANNIHILATION 

If DM is made up of a so-called "thermal relic" parti- 
clc""^, its abundance today is set by its annihilation cross 
section in the early universe. The thermal relic abun- 
dance calculation relating today's abundance of DM to 
the properties of the DM particle (its mass and anni- 
hilation cross section) is straightforward and elegant, 
and has been described in pedagogical detail previously 
^0, 61, 62] . We briefly summarize the story here. 

At sufficiently early times, the DM particles are in ther- 
mal equilibrium with the rest of the universe. As long as 
they remain relativistic (T 3> ti^), their creation and de- 
struction rates are balanced, and hence their co-moving 
abundance remains constant. Once the universe cools 
below the DM particle's rest-mass (T < m^), its equi- 
librium abundance is suppressed by a Boltzmann factor 
exp(— m^^/r). If equilibrium had been maintained until 
today, the DM particles would have completely annihi- 
lated away. Instead the expansion of the universe comes 
to the rescue and causes the DM particles to fall out of 
equilibrium once the expansion rate (given by II{a), the 



^ An alternative DM candidate is the axion, a non-thermal relic 
particle motivated as a solution to the strong CP problem j59j. 
Since it doesn't produce an annihilation signal today, we don't 
further consider it here. 



y ~ m^/20, with only a weak additional 
logarithmic dependence on the mass and cross section 
of the DM particle. A back of the envelope calculation 
results in the following relation between fi^, the relic 
mass density in units of the critical density of the uni- 
verse pcrit = 3ifo /SttG, and {av) , the thermally averaged 
velocity-weighted annihilation cross section: 



3 X 10- 



{(TV 



(1) 



Note that this relation is independent of m^. The 
WMAP satellite's measurement of the DM density is 
uj^ = 0.1131 ± 0.0034 [13, implying a value of 



3 X 10"^^ cm^ s-\ 



(2) 



A more accurate determination of (av) must rely on a 
numerical solution of the Boltzmann equation in an ex- 
panding universe, taking into account the full temper- 
ature dependence of the annihilation rate, including the 
possibilities of resonances and co-annihilations into other, 
nearly degenerate dark sector particles [e.g. [H, [6^1 ■ It is 
a remarkable coincidence that this value of (av) is close 
to what one expects for a cross section set by the weak 
interaction. This is the so-called "WIMP miracle" , and it 
is the main motivation for considering weakly interacting 
massive particles (WIMPs) as prime DM candidates. 

The Standard Model of particle physics actually pro- 
vides one class of WIMPs, massive neutrinos. Although 
neutrinos thus constitute a form of DM, they cannot 
make up the bulk of it, since their small mass, ^ < 
0.63 eV [il], implies a cosmological mass density of only 
cjjy = 7.1 X 10"'^. The attention thus turns to extensions 
of the Standard Model, which themselves are theoreti- 
cally motivated by the hierarchy problem (the enormous 
disparity between the weak and Planck scales) and the 
quest for a unification of gravity and quantum mechanics. 
The most widely studied class of such models consists of 
supersymmetric extensions of the Standard Model, al- 
though models with extra dimensions have received a 
lot of attention in recent years as well. Both of these 
approaches offer good DM particle candidates: the light- 
est supersymmetric particle (LSP), typically a neutralino 
in R-parity conserving supersymmetry, and the lightest 
Kaluza-Klein particle (LKP), typically the B^^^ parti- 
cle, the first Kaluza-Klein excitation of the hypercharge 
gauge boson, in Universal Extra Dimension models. For 
much more information, we recommend the comprehen- 
sive recent review of particle DM candidates by Bertone, 
Hooper & Silk 

The direct products of the annihilation of two 
DM particles are strongly model dependent. Typi- 
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FIG. 1: A schematic of the different sources and energy distri- 
butions of 7-rays from WIMP annihilation. Top: Secondary 
photons arising from the decay of neutral pions produced in 
the hadronization of primary annihilation products. Middle: 
Internal bremsstrahlung photons associated with charged an- 
nihilation products, either in the form of final state radiation 
(FSR) from external legs or as virtual internal bremsstrahlung 
(VIB) from the exchange of virtual charged particles. Bot- 
tom: Mono-chromatic line signals from the prompt annihila- 
tion into two photons or a photon and Z boson. This process 
occurs only at loop level, and hence is typically strongly sup- 
pressed. 



cal channels include annihilations into charged leptons 
(e+e", r+T^), quark- antiquark pairs, and gauge 
and Higgs bosons {W'^W' ,Z,h). In the end, how- 
ever, the decay and hadronization of these annihilation 
products results in only three types of emissions: (i) 
high energy neutrinos, (ii) relativistic electrons and pro- 
tons and their anti-particles, and (iii) 7-ray photons. 
Additional lower energy photons can result from the 
interaction of the relativistic electrons with magnetic 
fields (synchrotron radiation), with interstellar material 
(bremsstrahlung), and with the CMB and stellar radia- 
tion fields (inverse Compton scattering). In the follow- 
ing we will focus on the 7-rays, since they are likely the 
strongest signal from Galactic DM substructure. 7-rays 



are produced in DM annihilations in three ways (see ac- 
companying Fig. [1]) 

(i) Since the DM particle is neutral, there is no di- 
rect coupling to photons. Nevertheless, copious 
amounts of secondary 7-ray photons can be pro- 
duced through the decay of neutral pions, tt*^ — > 77, 
arising in the hadronization of the primary annihi- 
lation products. Since the DM particles are non- 
relativistic, their annihilation results in a pair of 
mono-energetic particles with energy equal to m^, 
which fragment and decay into 7r-meson dominated 
"jets" . In this way a single DM annihilation event 
can produce several tens of 7-ray photons. The re- 
sult is a broad spectrum with a cutoff around m^. 

(ii) An important additional contribution at high 
energies (E ^ to^) arises from the internal 
bremsstrahlung process [6^ , which may occur with 
any charged annihilation product. One can dis- 
tinguish between final state radiation, in which 
the photon is radiated from an external leg, and 
virtual internal bremsstrahlung, arising from the 
exchange of a charged virtual particle. Note 
that neither of these processes requires an exter- 
nal electromagnetic field (hence the name internal 
bremsstrahlung). The resulting 7-ray spectrum is 
peaked towards E ^ rriy. and exhibits a sharp cut- 
off. Although it is suppressed by one factor of the 
coupling a compared to pion decays, it can pro- 
duce a distinctive spectral feature at high energies. 
This could aide the confirmation of a DM annihila- 
tion nature of any source and might allow a direct 
determination of m^. 

(iii) Lastly, it is possible for DM particles to directly 
produce 7-ray photons, but one has to go to loop- 
level to find contributing Feynman diagrams, and 
hence this flux is typically strongly suppressed by 
two factors of a (although exceptions exist ^67]). 
On the other hand, the resulting photons would 
be mono-chromatic, and a detection of such a line 
signal would provide strong evidence of a DM an- 
nihilation origin of any signal. While annihilations 
directly into two photons, XX ~^ 11^ would pro- 
duce a narrow line a,t E = m^, in some models it 
is also possible to annihilate into a photon and a Z 
boson, XX ~* 7-^1 ^ind this process would result in 
a somewhat broadened (due to the mass of the Z) 
hne at E ^ "^xl^ ~ 

The relative importance of these three 7-ray produc- 
tion channels and the resulting spectrum dN^/dE de- 
pend on the details of the DM particle model under con- 
sideration. For any given model, realistic 7-ray spec- 
tra can be calculated using sophisticated and publicly 
available computer programs, such as the PYTHIA Monte- 
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Carlo event generator [68| , which is also contained in the 
popular DarkSUSY package [6^ . 



DARK MATTER SUBSTRUCTURE AS 
DISCRETE 7-RAY SOURCES 

DM subhalos as individual discrete 7-ray sources hold 
great potential for provi ding a "smo king gun" signature 



of DM annihilation [23, S 0, S H O, S H E , 
771 . . Compared to diffuse 7-ray annihilation 
signals, these discrete sources should be easier to distin- 
guish from astrophysical backgrounds and foregrounds 
[79] , since a) typical astrophysical sources of high energy 
7-rays, such as pulsars and supernova remnants, are very 
rare in dwarf galaxies, owing to their predominantly old 
stellar populations, b) the DM annihilation flux should 
be time-independent, c) angularly extended, and d) not 
exhibit any (or only very weak) low energy emission due 
to the absence of strong magnetic fields or stellar radia- 
tion fields. 

We can distinguish between DM subhalos hosting a 
Milky Way dwarf satellite galaxy and dark clumps that, 
for whatever reason, don't host a luminous stellar popu- 
lation, or one that is too faint to have been detected up 
to now. Before we go on to discuss the prospects of de- 
tecting a DM annihilation signal from these two classes 
of sources, we review the basic properties of DM subhalos 
common to both. 

Numerical simulations have shown that pure DM 
(sub)halos have density profiles that are well described 
by a Navarro, Frenk fc White (NFW) [sOl profile over a 
wide range of masses 
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PNFw{r) 



(r/r,)(l +r/r,)2 



(3) 



The parameter indicates the radius at which the log- 
arithmic slope 7(r) = = —2, and p{rs) = ps- The 
very highest resolution simulations have recently pro- 
vided some indications of a flattening of the density pro- 
file in the innermost regions 15, 83]. In this case a so- 
called Einasto profile may provide a better overall fit. 



PEinasto(?') = Ps exp 



a V V f . 



(4) 



Here the additional parameter a governs how fast the 
profile rolls over, and has been found to have a value 
of a « 0.17 ± 0.03 in numerical simulations [s^. Note 
that the two density profiles actually do not differ very 
much until r <C (cf . top panel of Figure [5]) . Sim- 
ulated DM halos are of course not perfectly spherically 
symmetric, and instead typically exhibit prolate or tri- 
axial iso-density contours that become more elongated 
towards the center [sl]. The degree of prolateness de- 
creases with mass, and galactic subhalos have axis ratios 
of > 0.7 fH. 
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FIG. 2; A comparison of NFW and Einasto (q = 0.17) radial 
profiles of density (top, dark lines, left axis), circular velocity 
(top, light lines, right axis), enclosed annihilation luminosity 
(bottom, dark lines, left axis), enclosed mass (bottom, light 
lines, right axis). The density profiles have been normalized 
to have the same Vmax and rvmax. 



The "virial" radius i?vir of a halo is defined as the 
radius enclosing a mean density equal to AvirPo, where 
Avir ~ 389 [86] and po is the mean density of the universe. 
The corresponding virial mass Mvir is the mass within 
i?vir, and a halo's concentration can then be defined as 
c = Rvir/rs- While these quantities are well defined for 
isolated halos and commonly used in analytic models, 
they are somewhat less applicable to galactic subhalos, 
since the outer radius of a subhalo is set by tidal trun- 
cation, which depends on the subhalo's location within 
its host halo. Furthermore, in numerical simulations it 
is difficult to resolve in low mass subhalos. For this 
reason we prefer to work with I^nax, the maximum of the 
circular velocity curve Vc(r)'^ = GM{< r)/r and a proxy 
for a subhalo's mass, and rvmax, the radius at which 
Knax occurs. Thcsc quantities are much more robustly 
determined for subhalos in numerical simulations than 
(M, c). Note that even (V^ax, ?'Vmax) can be affected by 
tidal interactions with the host halo, especially for sub- 
halos close to the host halo center. For this reason we 
also sometimes consider Vpcak, the largest value of Knax 
that a subhalo ever acquired during its lifetime (i.e. be- 
fore tidal stripping began to lower its I^nax) and rypcak, 
the corresponding radius. 

Since DM annihilation is a two body process, its rate 
is proportional to the square of the local density, and the 
annihilation "luminosity" is given by the volume integral 
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of p(r)2, 



£(< r) 



dV. 



(5) 



C has dimensions of (mass)^ (length)"^, and we express 
it in units of Mq pc~^. In order to convert to a conven- 
tional luminosity, one must multiply by a particle physics 
term, 



L = 



2 {<^V) 



(6) 



where is the mass of the DM particle and (aw) the 
thermally averaged velocity-weighted annihilation cross 
section discussed in the previous section. This is the 
total luminosity, but we are interested here only in the 
fraction emitted as 7-rays. Furthermore, a given detector 
is only sensitive to 7-rays above a threshold energy of i?th 
and below a maximum energy of -Emax- In that case the 
effective 7-ray luminosity is 



cff 



2m2 



E- 



Ie 



'-dE 



(7) 



where dN^/dE is the spectrum of 7-ray photons pro- 
duced in a single annihilation event. 

A comparison of the enclosed luminosity and mass pro- 
files is shown in the bottom panel of Figure O Clearly, C 
is much more centrally concentrated than M : ^ 90% 
of the total luminosity is produced within r^, com- 
pared with only 10% of the total mass. In terms of 
(14nax7 ^Vmax)j the total luminosity of a halo is given by 



1/4 



(8) 



where / is an 0(1) numerical factor that depends on the 
shape of the density profile; for an NFW profile / = 
1.227, and for an a = 0.17 Einasto profile / = 1.735. In 
physical units, the total annihilation luminosity is 



£ = 1'^ X 10^ M| pc- 



20 km s- 



Ikpc 



(9) 

for NFW (top) and a = 0.17 Einasto (bottom). Note 
that even though the slope of the Einasto profile is 
shallower than NFW in the very center, the total lu- 
minosity exceeds that of an NFW halo with the same 
(Vinax, ?'Vmax)- This is duc to the fact that the Einasto 
profile rolls over less rapidly than the NFW profile, and 
actually has slightly higher density than NFW between 
Ts and a cross-over point at ^ lO^^rg. 



MILKY WAY DWARF SPHEROIDAL GALAXIES 

There are several advantages of known dwarf satellite 
galaxies as DM annihilation sources: firstly, the kinemat- 
ics of individual stars imply mass-to-light ratios of up to 



several hundred 87, and hence there is an a 

priori expectation of high DM densities; secondly, since 
we know their location in the sky, it is possible to directly 
target them with sensitive atmospheric Cerenkov tele- 
scopes (ACT) such as H.E.S.S., VERITAS, MAGIC, and 
STACEE, whose small field of view makes blind searches 
impractical; lastly, our approximate knowledge of the dis- 
tances to many dwarf satellites would allow a determina- 
tion of the absolute annihilation rate, which may lead 
to a direct constraint on the annihilation cross section, 
if the DM particle mass can be independently measured 
(from the shape of the spectrum, for example). 

Recent observational progress utilizing the Sloan Digi- 
tal Sky Survey (SDSS) has more than doubled the num- 
ber of known dwarf spheroidal (dSph) satellite galaxies 
orbiting the Milky Way 0, i, i, 0, i, 0, 0, [III, [ll , rais- 
ing the total from the 9 "classical" ones to 23. Many of 
the newly discovered satellites are so-called "ultra-faint" 
dSph's, with luminosities as low as 1,OOOL0 and only 
tens to hundreds of spectroscopically confirmed mem- 
ber stars. Simply accounting for the SDSS sky coverage 
(about 20%), the total number of luminous Milky Way 
satellites can be estimated to be at least 70. Taking into 
account the SDSS detection limits (lH and a radial dis- 
tribution of DM subhalos motivated by numerical sim- 
ulations, this estimate c an g row to several hundreds of 
satellite galaxies in total 9^, l93| . 

In order to assess the strength of the DM annihila- 
tion signal from these dSph's, it is necessary to have an 
estimate of the total dynamical mass, or at least 14iax, 
of the DM halo hosting the galaxies. Owing to the ex- 
treme faintness of these objects and their lack of a de- 
tectable gaseous component [s^, it has been very dif- 
ficult to obtain kinematic information that allows for 
such measurements. Progress has been made through 
spectroscopic observations of individual member stars, 
whose line-of-sight velocity dispersions have confirmed 
that these objects are in fact strongly DM dominated 
[stL IssL [sil . [90 1 . Such data best constrain the enclosed 
dynamical mass within the stellar extent, which on av- 
erage is about 0.3 kpc for current data sets. A recent 
analysis has determined A/0.3 = M{< 0.3 kpc) for 18 of 
the Milky Way dSph's, and found that, surprisingly, they 
all have M0.3 ~ 10^ Mq to within a factor of two [3l| . 

State-of-the-art cosmological numerical simulations of 
the formation of the DM halo of a Milky Way scale 
galaxy, such as those of the Via Lactea Project 13, T^l 
and the Aquarius Project [lj| , have now reached an ad- 
equate mass and force resolution to directly determine 
M0.3 in their simulated subhalos. This makes it possible 
to infer the most likely values of (Kiax, ?'Vmax) for a Milky 
Way dSph of a given M0.3 and Galacto-centric distance 
by identifying all simulated subhalos with comparable 
Mo. 3 and D and averaging over their (V^ax, fvmax)- This 
analysis was performed for the 9 "classical" dwarfs using 
the Via Lactea I simulation 22!1 , and we extend it here to 
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TABLE I: The pro perties of likely DM subhalos of the 18 Milky Way dSph galaxies for which M0.3 values (column 3) have 
been published [33]. Vmax and rvmax are the maximum circular velocity and its radius, Vpeak and rvpeak the largest Knax a 
subhalo ever acquired and its corresponding radius. The first number is the median over all Via Lactea II subhalos matching 
the dSph's distance and Mo. 3, the numbers in parentheses the 16"^ and 84**^ percentiles. (See text for details.) 



all 18 dwarfs published in j3l| and with the more recent 
and higher resolution Via Lactea II (VL2) simulation. 

We randomly generated 100 observer locations at 8 
kpc from the VL2 host halo center, and selected, for 
each Milky Way dSph in [lH separately, all simulated 
subhalos with distances within 40% and numerically de- 
termined M0.3 within the published 1 — a error bars. We 
then determined the median value and the 16**^ and 84'^ 
percentiles of (Knax, fvmax) and (Vpcak, ^vpcak) for each 
dSph. These values are given in Table HI The median 
values of Knax range from 5.0 km s~^ (Leo IV) to 19 km 
s~^ (Draco, Leo I), and of Fpoak from 6.7 km s~^ (Leo 
IV) to 30 km s^^ (Ursa Minor, Carina). Note that, as 
expected, dSph's closer to the Galactic Center typically 
show a larger reduction from Ipcak to Vmax, sometimes 
by more than a factor of 2. 

In the same fashion, we then determine the most 
likely annihilation luminosities for the 18 dSph's by using 
Eq. ([9]) for an NFW profile to calculate the total lumi- 
nosity £tot for every simulated subhalo. Additionally we 
also determine £0.3 j the luminosity within 0.3 kpc from 
the center, motivated by the fact that we only have dy- 
namical evidence for a DM dominated potential out to 



this radius. Lastly we also consider two measures of the 
brightness of each halo: JFtot — 'Ctot/47rZ3^, the total ex- 
pected flux from the dSph, and J^c, the flux from a cen- 
tral region subtending 0.15°, which is comparable to the 
angular resolution of Fermi above 3 GeV. J-c thus corre- 
sponds to the brightest "pixel" in a Fermi 7-ray image 
of a subhalo. These numbers are given in Table [TTl 

Current observational constraints 

Several ACT have performed observations of a handful 
of dSph's. 

• The H.E.S.S. array (consisting of four 107 m^ tele- 
scopes with a 5° fleld of view and an energy thresh- 
old of 160 GeV [il]) has obtained an 11 hour ex- 
posure of the Sagittarius dwarf galaxy. No 7-ray 
signal was detected, resulting in a flux limit of 
3.6 X 10-^2 cm-2 s-i (95% confldence) at E > 250 
GeV, and a corresponding limit on the cross section 
of (av) < 10^^^ cm^ s^^ for an NFW profile and 
(cw) ^ 2 X 10"^'°^ cm'^ s~^ for a cored profile (for a 
= 100 GeV - 1 TeV neutrahno) fii]. Note that 
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TABLE II: Estimated luminosities and fluxes for the 18 dSph from Table U] £tot is the total luminosity and £0.3 the luminosity 
from within the central 0.3 kpc. Ttot = £tot/47r_D^ is the total flux and Ta the flux from a central region subtending 0.15° 
(about the angular resolution of Fermi above 3 GeV). The first number is the median over all subhalos matching the dSph 
distance and Mo. 3, the numbers in parentheses are the 16**^ and 84*'' percentiles. The table is ordered by decreasing Ttot. 



the Sagittarius dwarf is undergoing a strong tidal 
interaction with the Milky Way galaxy [96] , and no 
confident determination of A/0.3 has been possible. 

• The VERITAS array (consisting of four 144 
telescopes with a 3.5° field of view and an energy 
threshold of 100 GeV [97l|) has conducted a 15 hours 
observation of Willman 1 and 20 hour observations 
each of Draco and Ursa Minor [45|. No 7-ray sig- 
nal was detected at a flux limit of 1% of the flux 
from the Crab Nebula, corresponding to a limit of 
2.4 X 10-12 g-i ^95^^ confldence) at > 200 

GeV fal. 



• Additionally the MAGIC and STACEE fi^ 

telescopes have reported observations of Willman 1 
and Draco, resulting in comparable or slightly 
higher flux limits. 

To convert the values of Tc in TablelTTlinto physical fluxes 
that can directly be compared to these observational lim- 
its, it would be necessary to obtain values of the particle 
physics term of Eq. ([7]) by performing a scan of the DM 
model parameter space. This is beyond the scope of this 



work, but a si milar analy sis has been performed by oth- 
ers [nl, [ii m EoS Eoif. current ACT observations of 
dSph's are beginning to directly constrain DM models, 
and future longer exposure time observations of addi- 
tional dSph's (in particular Segue 1 and Ursa Major II) 
with a lower threshold energy hold great potential. We 
also eagerly await the first Fermi data on fluxes from the 
known dSph galaxies. 



DARK CLUMPS 

An annihilation signal from dark clumps not associ- 
ated with any known luminous stellar counterpart would 
provide evidence for one of the fundamental implications 
of the CDM paradigm of structure formation: abundant 
Galactic substructure. Barring a serendipitous discovery 
with an ACT, the discovery of such a source will have to 
rely on all-sky surveys, such as provided by Fermi. Of 
course even a weak and tentative identification of a dark 
clump with Fermi could be followed up with an ACT. 

Unlike for known dSph galaxies, for which we at least 
have some astronomical observations to guide us, we 
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FIG. 3: The annihilation flux .T^tot from subhalos in the Via Lactea II simulation versus their Mo. 3, Vmax, and Vpcak- The 
gray shaded areas indicate regions containing subhalos with Ttot as least as high as the fifth-brightest Milky Way dSph galaxy 
(Carina), but with Mo. 3, Knax, Vpeak below that of the known dSph's, i.e. probable dark clumps. Only one of the 100 random 
observer locations used in the analysis is shown here. 



must rely entirely upon numerical simulations to quantify 
the prospects of detecting the annihilatio n si gnal from 
dark clumps. Recent significant progress 102l | notwith- 
standing, it is at present not yet possible to perform real- 
istic cosmological hydrodynamic galaxy-formation simu- 
lations, which include, in addition to the DM dynamics, 
all the relevant baryonic physics of gas cooling, star for- 
mation, supernova and AGN feedback, etc. that may 
have a significant impact on the DM distribution at the 
centers of massive halos. Instead we make use of the 
extremely high resolution, purely coUisionless DM-only 
Via Lactea II (VL2) simulation [l3|, which provides an 
exquisite view of the dumpiness of the Galactic DM dis- 
tribution, but at the expense of not capturing all the rel- 
evant physics at the baryon-dominated Galactic center. 
For the abundance, distribution, and internal properties 
of the DM subhalos that are the focus of this work, the 
neglect of baryonic physics is less of a problem, since they 
are too small to allow for much gas cooling and significant 
baryonic effects (this is supported by the high mass-to- 
light ratios observed in the Milky Way dSph's), although 
tidal interactions with the Galactic stellar and gaseous 
disk might significantly affect the population of nearby 
subhalos. 

With a particle mass of 4, 100 M© and a force softening 
of 40 pc, the VL2 simulation resolves over 50,000 subha- 
los today within the host's r2oo — 402 kpc (the radius 
enclosing an average density 200 times the mean matter 
value). Above ^ 200 particles per halo, the differen- 
tial subhalo mass function is well-fit by a single power 
law, dN/dM ~ M~^-^, and the cumulative Vmax func- 
tion is N{> Knax) ^ V-^^ The radial distribu- 
tion of subhalos is "anti-biased" with respect to the host 



halo's density profile, meaning that the mass distribu- 
tion becomes less clumpy as one approaches the host's 
center l^, 8^ • Similar results have been obtained by the 
Aquarius group [flll^l- Typical subhalo concentrations, 
defined as Ay = (p(< ''Vmax))/Pcrit, grow towards the 
center, ow ing to a combination of earlier formation times 
IO3I 104 1 and stronger tidal stripping of central subha- 



los: VL2 subhalos on average have a 60 times higher Ay 
at 8 kpc than at 400 kpc [l3|. Note that this also im- 
plies ~ 7 times higher annihilation luminosities for cen- 



tral subhalos, since £ ^ Knax/'^V: 



' m 



The 



counter-acting trends of decreasing relative abundance of 
subhalos and increasing annihilation luminosity towards 
the center makes it more difficult for (semi-)analytical 
methods to accurately assess the role of subhalos in the 
Galactic annihilation signal, and motivate future, even 
higher resolution, numerical simulations of the formation 
and evolution of Galactic DM (sub-) structure. A direct 
analysis of the VL2 simulations in terms of the detectabil- 
ity with Fermi of individual subhalos was performed by 
[78| . They found that for reasonable particle physics pa- 
rameters a handful of subhalos should be able to outshine 
the astrophysical backgrounds and would be detected at 
more than Scr significance over the lifetime of the Fermi 
mission. 

As discussed in the previous section, we have directly 
calculated the annihilation luminosities for all VL2 sub- 
halos using Eq. ([9]) and assuming an NFW density pro- 
file. The luminosities would be ~ 40% higher if an 
Einasto (a = 0.17) profile had been adopted instead. We 
then converted these luminosities to fluxes by dividing 
by 47rZ?^, where the distances D were determined for 100 
randomly chosen observer locations 8 kpc from the host 
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halo center. The resulting values of Ttot are plotted in 
Figure [3l for just one of the 100 observer positions, as 
a function of the subhalos' M0.3, Knax, and 14>eak- Al- 
though the distributions show quite a bit of scatter, in 
all three cases a clear trend is apparent of more massive 
subhalos having higher fluxes. This trend could simply 
be the result of the higher luminosities of more massive 
halos, but one might have expected smaller mass subha- 
los to be brighter, since their greater abundance should 
result in lower typical distances and hence higher fluxes. 
This latter effect could be artificially suppressed in the 
numerical simulations, if smaller mass subhalos, whose 
dense centers are not as well resolved, were more eas- 
ily tidally disrupted closer to the Galactic Center, or if 
the subhalo finding algorithm had trouble identifying low 
mass halos in the high background density central region. 
In Figure |4] we plot the subhalos' Vmax against their dis- 
tance to the host halo center D. There appears to be a 
dearth of the lowest Knax subhalos (Vmax ^ 2 km s~^) at 
small distances, but at the moment it is not clear whether 
this suppression is a real effect or a numerical artifact. 
It's also worth noting that such small subhalos might 
be more susceptible to disruption by interactions with 
the Milky Way's stellar and gaseous disk. At any rate, 
we can obtain an analytic estimate of the scaling of the 
typical subhalo flux with Knax by noting that the lumi- 
nosity scales as £ ~ V^^^^fKv and the typical distance 
as D ^ ~ y^/^ (since dn/dKiax - ^^nax)- The 

typical fiux should thus scale as ~ £/Z?^ ~ V^mxVAv, 
and would be higher for more massive subhalos at a fixed 
Av. Actually lower Vinax subhalos might be expected to 
have higher Ay due to their earlier formation times, but 
it remains to be seen to what degree this expectation is 
borne out in numerical simulations. 

The points in Figure [3] can be directly compared with 
the values for the known Milky Way dSph's in Tables H] 
and nil it appears that there are many DM subhalos at 
least as bright as the known Milky Way dSph's. This 
impression is confirmed by the top panel of Figure [5l in 
which we show the cumulative number of subhalos with 
fiuxes greater than ^tot and Tc- These distributions were 
obtained by averaging over 100 randomly chosen observer 
locations 8 kpc from the host halo center. The mean 
number of DM subhalos with Ttot greater than that of 
(Carina, Draco, Ursa Minor, Ursa Major, Segue 1) is (90, 
54, 43, 17, 13), and the corresponding numbers for J^c are 
(96, 62, 49, 24, 19). This demonstrates that if a DM anni- 
hilation signal from any of the known Milky Way dSph's 
is detected, then many more DM subhalos should be vis- 
ible. The plot also implies that Segue 1, the dSph with 
the highest J-tot and J^c of the currently known sample, 
is unlikely to be the brightest DM subhalo in the sky. 
Of course some of these additional bright sources could 
very well have stellar counterparts that have simply been 
missed so far, due to the limited sky coverage of current 
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D [kpc] 

FIG. 4: VL2 subhalo Vmax vs. distance to host halo D. 




FIG. 5: Top: The cumulative number of subhalos with flux 
exceeding Ttot, J~c- Bottom: The fraction of dark clumps, 
i.e. subhalos likely not hosting any stars and defined by 
M0.3 < 5 X 10*^ M0, Vmax < 8 km s~\ or Vpcak < 14 km s"\ 
as a function of Umiting flux J-toi- These distributions are 
averages over 100 randomly chosen observer locations 8 kpc 
from the host halo center. 
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surveys or insufficiently deep exposures. To assess what 
fraction of high flux sources are likely to be genuinely 
dark clumps without any stars, we split the sample by a 
limiting value of M0.3 = 5 x 10^ Mq, T^iax = 8 km s~^, 
and Knax = 14 km s~^. We assume that DM subhalos 
below these limits are too small to have been able to 
form any stars, and hence are truly dark clumps. Of the 
known dSph's listed in Table |T] only Leo IV falls below 
these limits. In the bottom panel of Figure [5] we plot 
/dark(> -^tot), the fractiou of subhalos without stars, as 
a function of the limiting annihilation flux J- tot- /dark 
falls monotonically with J^tot, which makes sense given 
that higher flux sources are typically more massive and 
hence more likely to host stars. Between 30 and 40% of 
all DM subhalos brighter than Carina are expected to be 
dark clumps. This fraction drops to 10% for subhalos 
brighter than Segue 1. 



Boost Factor? 

The analysis presented here so far has been limited to 
known dSph galaxies and clumps resolved in the VL2 
simulation, whose resolution limit is set by the avail- 
able computational resources, and has nothing to do with 
fundamental physics. Indeed, the CDM expectation is 
that the dumpiness should continue all the way down to 
the cut-off in the matter power spectrum, set by colli- 
sional da mping and free streaming in the early universe 
lolIiQi- For typical WIMP DM, this cut-off occurs at 



masses of mo = 10"^^ to 10"'' M© [25|,|26l, some 10 to 20 
orders of magnitude below VL2's mass resolution. Since 
the annihilation rate goes as and (p^) > (p)^, this sub- 
resolution dumpiness will lead to an enhancement of the 
total luminosity compared to the smooth mass distribu- 
tion in the simulation. 

The magnitude of this so-called substructure boost fac- 
tor depends sensitively on the properties of subhalos be- 
low the simulation's resolution limit, in particular on the 
behavior of the concentration-mass relation. A simple 
power law extrapolation of the contribution of simulated 
subhalos to the total luminosity of the host halo leads to 
boosts on order of a a few hundred [H^] . More sophisti- 
cated (semi-) analytical models, accounting for different 
possible continuations of the concentration-mass relation 
to lower masses, t ypically f ind smaller boosts of around 
a few tens ^ m}^, [iQI | . 



More importantly, this boost refers to the enhancement 
of the total annihilation luminosity of a subhalo, but this 
is not likely the quantity most relevant for detection. At 
the distances where subhalos might be detectable as in- 
dividual sources, their projected size exceeds the angular 
resolution of today's detectors. The surface brightness 
profile from annihilations in the smooth DM component 
would be strongly peaked towards the center (yet prob- 
ably still resolved by Fermi [zl]), owing to the p(r)^ de- 



pendence of the annihilation rate. The luminosity contri- 
bution from a subhalo's sub-substructure population (i.e. 
its boost), however, is much less centrally concentrated: 
at best it follows the subhalo's mass density profile p{r), 
although it might very well even be anti-biased. This 
implies that substructure would preferentially boost the 
outer regions of a subhalo, where the surface brightness 
typically remains below the level of astrophysical back- 
grounds and hence doesn't contribute much to the detec- 
tion significance. In other words, the boost factor might 
apply to J-tot, but much less (or not at all) to Tc] yet it 
is Tc that is likely to determine whether a given subhalo 
can be detected with Fermi or an ACT. It thus seems un- 
likely that the detectability of Galactic subhalos would 
be significantly enhanced by their own substructure^ . On 
the other hand, a substructure boost could be very im- 
portant for diffuse DM annihilation signals, either from 
extragalactic sources, where the boost would simply in- 
crease the overall amplitude |l08l |. or from Galactic DM, 
where the boost could affect the amplitude and angular 
profile of the signal, as well as the power spectrum of its 



anisotropics ^ 0, [Hll, [H, [zl, [ji 



SUMMARY AND CONCLUSIONS 

In this work we have reviewed the DM annihilation sig- 
nal from Galactic subhalos. After going over the basics 
of the annihilation process with a focus on the resulting 
7-ray output, we summarized the properties of DM sub- 
halos relevant for estimating their annihilation luminos- 
ity. In the remainder of the paper we used the Via Lactea 
II simulation to assess the strength of the annihilation 
flux from both known Galactic dSph galaxies as well as 
from dark clumps not hosting any stars. By matching the 
distances D and central masses M0.3 of simulated subha- 
los to the corresponding published values of 18 known 
dSph's, we were able to infer most probable values, and 
the 1-a scatter around them, for Vmax and rvmax, and 
hence for the annihilation luminosity C and fiux of all 
dwarfs. According to this analysis, the recently discov- 
ered dSph Segue 1 should be the brightest of the known 
dSph's, followed by Ursa Major II, Ursa Minor, Draco, 
and Carina. Further, we showed that if any of the known 
Galactic dSph's are bright enough to be detected, then 
at least 10 times more subhalos should appear as visi- 
ble sources. Some of these would be as-of-yet undiscov- 
ered luminous dwarf galaxies, but a significant fraction 
should correspond to dark clumps not hosting any stars. 



^ This is in contrast to many previous claims in the literature, in- 
cluding some by the present author [e.g. [T^. A re-analysis of 
that work (in progress) with an improved treatment of the angu- 
lar dependence of the substructure boost, indeed finds that the 
boost only weakly increases the number of detectable subhalos. 
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The fraction of dark clump sources is 10% for subhalos 
at least as bright as Segue 1 and grows to 40% for sub- 
halos brighter than Carina. Lastly, we briefly considered 
the role that a substructure boost factor should play in 
the detectability of individual Galactic dSph's and other 
DM subhalos. We argued that any boost is unlikely to 
strongly increase their prospects for detection, since its 
shallower angular dependence would preferentially boost 
the outer regions of subhalos, which typically don't con- 
tribute much to the detection significance. 

Several caveats to these findings are in order. Prob- 
ably the most important of these is that our simulation 
completely neglects the effects of baryons. Gas cooling, 
star formation, and the associated feedback processes are 
unlikely to strongly affect most subhalos, owing to their 
low masses. However, tidal interactions with the bary- 
onic components of the Milky Way galaxy might do so. 
The Sagittarius dSph, for example, is thought to be in 
the process of complete disruption from tidal interactions 
with the Milky Way. A second caveat is that our analysis 
is based on only one, albeit very high resolution, numer- 
ical simulation, and so we cannot assess the importance 
of cosmic variance, or the dependence on cosmological 
parameters such as erg and n^. Oth er work h as found 
considerable halo-to-halo scatter 14, llOOl . IllOj . with a 
factor of ~ 2 variance in the total subhalo abundance, 
for example. 

These caveats motivate further study and future, 
higher resolution numerical simulations, including the ef- 
fects of baryonic physics. The characterization of the 
Galactic DM annihilation signal is of crucial importance 
in guiding observational efforts to shed light on the na- 
ture of DM. We are hopeful that in the next few years the 
promise of a DM annihilation signal will come to fruition, 
and will help us to unravel this puzzle. 
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